Single-electron transport through the vortex core levels in clean superconductors 
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We develop a microscopic theory of single-electron transport in N-S-N hybrid structures in the 
presence of applied magnetic field introducing vortex lines in a superconductor layer. We show 
that vortex cores in a thick and clean superconducting layer are similar to mesoscopic conducting 
channels where the bound core states play the role of transverse modes. The transport through not 
very thick layers is governed by another mechanism, namely by resonance tunneling via vortex core 
levels. We apply our method to calculation of the thermal conductance along the magnetic field. 
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Electron transport through various hybrid structures 
is in the focus of current nanoscale physics research. Of 
a special interest are normal metal (N) - superconductor 
(S) - normal metal (N) trilayers where a superconduct- 
ing gap Ag suppresses single-particle transport, making 
charge transfer transparency very sensitive to the exter- 
nal controlling parameters. The electron transmission 
through such N-S-N structure at low energies is associ- 
ated with two-particle Andreev processes. If the thick- 
ness d of the superconducting slab is much larger than the 
coherence length ^ the electrons incident on the slab are 
reflected as holes, and the normal current converts into 
the supercurrent. Single-electron tunneling through an 
N-S-N structure decays exponentially with the slab thick- 
ness giving rise, in particular, to the exponential drop off 
of the electronic contribution to the thermal conductance 
P, • A single-particle transport through NSN recovers 
by applying a strong magnetic field, which creates vortex 
lines where the gap in the spectrum is suppressed. In the 
present Letter we develop a regular theoretical descrip- 
tion of the low-energy single-particle transfer through the 
vortex core states in clean type II superconductors. Since 
the single particle contribution to electric conductivity is 
short-circuited by supercurrent, we focus on the thermal 
conductivity which in this case is the experimentally ac- 
cessible characteristic of the one-electron transport. 

Taking the simplest view of a vortex core as a normal 
channel with a density of states as in the normal metal, 
we arrive at a single-electron Sharvin conductance of a 
normal wire with the radius of the order of ^: 



where kp = pp/fi is the Fermi wave vector. A little bit 
more attentive second thought shows that in clean super- 
conductors the only trajectories that do not hit vortex 
core "walls" contribute to a single-particle conductivity. 
Indeed, if a vortex core is a normal cylinder surrounded 
by a superconductor, an electron flying into the core 
boundary is Andreev reflected as a hole back along its 



incidental path. Thus the single-electron transport along 
such a trajectory is blocked, and the only contribution to 
the conductance of the "Andreev wire" comes from tra- 
jectories that traverse freely the normal region. Within 
this model, such trajectories must have incident angles 
9 < ^/d, i.e., are confined to a solid angle ^^/<P. This 
would result in G^e oc {e^ /TTh){kFO^i^/df . Consid- 
ering below the realistic gap profile inside the core in the 
framework of the developed approach we will show that 
the drop of conductance with the thickness d is even more 
rapid and proportional to d~^. 

However, the low energy transport associated with the 
Caroli-deGennes-Matricon (CdGM) states 4^ propagat- 
ing along the core should obviously saturate for very thick 
superconducting slabs. The CdGM spectrum e^i[kz) as a 
function of the quantized angular momentum = n+1/2 
varies from — Aq to -l-Ag crossing zero as the impact pa- 
rameter b = —fi/kj. of the particle varies from —oo to 
+00. Here fc^ = — k1 is the wave vector in the 

plane perpendicular to the vortex, (r, (p, z) is a cylin- 
drical coordinate system with the z- axis chosen along 
the vortex line. For small /i the spectrum is t^{kz) — 
— /zAo/(fcr^). Transport carried by the quantized trans- 
verse modes is described by the Landauer formula. In 
the limit ^^/{kpS) ^ T ^ Aq, the number of modes is 
~ (r/Ao)fcFCj Eind for the resonant contribution to the 
single-particle conductance one gets 0: 



(2) 



(1) Here /i numerates the transverse modes with transparen- 
cies Tn open in the vortex core. This estimate can also 



be obtained from Eq. provided the group velocity is 
taken as Vg = hT^de^/ dkz ~ €^/hkF instead of the veloc- 
ity ~ I'F as in a normal tube. The conductance of Eq. 
(0) is by a factor (A:fO~^ ^ 1 smaller than that of Eq. 

One thus expects that the single-particle Andreev- 
wire conductance of the vortex core would transform into 
Eq. with increasing d. 
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As we already mentioned, the single-electron transport 
determines the behavior of thermal conductivity. The 
Wiedemann-Franz law would result in k cx TG^e/e^ for 
the thermal conductance along the vortex line. In the 
vicinity of Hc.2 , the thermal conductivity has been stud- 
ied theoretically in a number of papers (see, for instance, 
[M IE Q)' In dirty superconductors, the electron con- 
tribution to the thermal conductance along the vortices 
for a small concentration of vortex lines ja] agrees con- 
ceptually with Eq. QJ: n{B) ~ {B / Hc2)kn , where kn 
is the electron thermal conductivity in the normal state. 
Unfortunately, in clean superconductors this siniple esti- 
mate fails to describe the experimental data |9|, [lOj : the 
thermal conductance in the magnetic field direction ap- 
pears to be two orders of magnitude smaller. It was noted 
first in that this obvious conflict can be caused by 
a very small group velocity of the CdGM modes as dis- 
cussed above. Analysis of quantum transport through 
individual vortices is of particular importance for under- 
standing the properties of mesoscopic superconductors. 
The exotic vortex states in these systems are nowadays 
the focus of a considerable attention [ill IT3 . Il3| . 

Hereafter we concentrate on the low-energy single- 
particle transport through vortex cores in clean {£ d) 
type II superconductors in the low-field limit of sepa- 
rated vortices and develop a systematic approach for cal- 
culation of the thermal conductance along vortex cores. 
We study the transmission of an electron wave incident 
on the superconducting slab placed between two bulk 
normal metal electrodes assuming ideally transparent 
boundaries and neglecting the normal scattering. Con- 
sidering two extremes of infinite and finite slab thick- 
nesses we confirm the intuitive picture discussed above: 
For a not very thick slab, the transmission is determined 
by the semi-classical resonant tunneling through the en- 
ergy gapped region when the energy coincides exactly 
with one of the levels in the vortex core. The transmis- 
sion is proportional to the large Sharvin conductance; 
however, it decays exponentially with the slab thick- 
ness except for the trajectories that go almost parallel to 
the vortex axis. The exponential decay is thus replaced 
with a i^/d)^ power law, which gradually transforms 
the Andreev-like thermal conductance into a thickness- 
independent expression Eq. ((JJ as c? increases. 

Wave functions - Quasiparticles in the superconductor 
obey the Bogolubov-de Gennes (BdG) equations: 
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{-ihV + -A] - Ef 
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+ Av = eu , (3) 
V - A*u = -ev . (4) 



The wave vector kz along the vortex axis is a good quan- 



tum number, u — e '^Uk^, v — '^Vk^- We put 
A ^ |A|e*'^, and Uk, = e''^/'^+''"l'U , Vk, = e-'-^/^-i-v^]/. 
Equations (|3J|, l@J have eigenvalues e^(fcz) for the CdGM 



bound states. We establish first a simple identity for the 
localized states. Calculating the derivative with respect 
to kz from the both sides of Eqs. ©, Q and using the 
normalization of the CdGM wave functions we find 



Ufj.k, 



hk^. 



'-A, 



v^k. 



d^r ■ 



h dkz 



(5) 



There is a considerable cancellation in Eq. (jSJ: the r.h.s. 
is by a factor of (A;f'C)~^ smaller than each term in the 
l.h.s. This is why Eq. ^ gives much smaller conductance 
than Eq. We will use Eq. lO later to derive the 

thermal conductance for a thick slab. 

However, the cancellation does not take place for a 
finite-thickness slab where the CdGM states are not truly 
localized. To consider this in more detail we use a semi- 
classical approach and look for U — {U, V) in the form 



= H^^^ [Kr) + {krr) w 



;,(-) 



r(i,2) 
assuming that w 



where H^'^'^' are the Hankel functions, / = 



1/4, 

(wi, W2) are slow functions of r and 

6^, 6 = —ii/kr and 



z. Following |l4l| we put x — \fr^ 
define the trajectories dx = ±(is± sin 0, dz — ds± cos 9 
for respectively, where k^ — kFCOsO and ds± is 

the distance along the corresponding trajectory. For a 
point {x, z) on the trajectory z = zq ± x cot 6 we obtain 



ih^kr dw^^'> 
m dx 



h^krb 



2m(x2 -I- 62) 
+iay\A\w^^'^ 







(6) 



For the functions z«^^^(x, z), the limit x ^ 00 corre- 
sponds to a wave radiating from the vortex into the 



bulk while 



„(-) 



{x, z) corresponds to a wave incident 
The condition of regularity at r = 



on the vortex. 

requires (0,z) = w^^^(0, z) at the classical turn- 
ing point, a; = 0. Let us put w''''-' (x, z) ~ w{x,z), 
w'^)(a;,z) — w{—x,z). The functions wi^2{x,z) satisfy 
Eq. ijBJ with the upper sign along the entire x axis. 
We introduce new functions ij and C through wi ~ 



W2 

drj 
dx 

dc 

dx 



gC "?/2 Qj^^ arrive at the equations 



2me 



2m|A| 



(a;2 + 62) n-^kr 



■ cos 77 



-m|A|/i ^fc^ ^ sinr/ 



(7) 
(8) 



The requirement that w vanishes at x ±00 is 77 = 
±7r/2 — (e/|Ao) + 2Trk. These values, however, are not 
stable for a general choice of the integration constants. 
A general solution for not very large positive x is 



77 



arctan ■ 



2m 



2K{x) 



L 



3 




FIG. 1: The coordinate dependence of 77 for 7 + ?7o = 0. The 
full line is for 7 > while the dashed line is for 7 < 0. 



where 770 is a constant and 

K{x) = mh-^k;^ I \A{x')\dx' 



For ?7o <C 1 and e ^ Aq, the function 77 is close to 7r/2 
for b ^ X < xq where xq ~ ^ln(l/|7|) and 
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2'mh 



[e-{b/x)\A\]e-^^'^^Ux . (10) 



7 measures the distance from a CdGM level: 7 = ex- 
actly when e = e^{kz). The validity of Eq. ^ is re- 
stricted by the condition 7^1 which generally holds 
if e ^ Aq. The function rj grows with |x| at distances 
X > Xq. To find its behavior in the region 77 ~ 1 we 
can neglect small terms with e and b/x in Eq. Q. The 
solution is 



tan 



(V _ ^\ 
V2 4/ 



Ce 



2K{x) 



(11) 



Matching with Eq. © at ^ ^ x ^ xq gives C — 
in ~^Vo) /2- For 7 -I- ?7o > 0, the function ry — ^ 37r/2 
while w diverges exponentially as x 00. If 7 + r/o < 0, 
the function 77 approaches — 7r/2, and w diverges again. 
However, if 7 -I- 770 = 0, the value 77 = 7r/2 is stable (see 
Fig. ^ and the wave function decays for x ^ 00. 

The solution at |x| <C xq for negative x is obtained from 
Eq. (O by replacing K{x) with —K{x). The function 77 
is close to — 7r/2 for 6 <C |x| ^ xq. Its behavior for 
\x\ > xo is determined by Eq. ifTT)) where C = 2/ (7 — 770). 
Equation then shows that — 37r/2 < rj{~oo) < —tt/2 
if 7 — 770 > 0. The function 77 grows and approaches 
r]{—oo) « — 37r/2 while w diverges. Similarly, —tt/2 < 
77(— 00) < n/2 if 7 — 770 < so that w also diverges. The 
value 77 — — 7r/2 is stable if only 7 — 770 = 0. The wave 
function thus decays at both ends if 7 = 770 = 0, which 
corresponds to a standard CdGM state. 

Reflection and transmission probabilities - For a super- 
conducting slab with a thickness d we thus have a linear 
combination of two solutions w = A^w'' + A^w^ where 



Ay and A< are constants. The first solution 2 h^^s 
^70 = —7 and decays at x ^ 00. Thus 77>(-t-|x|) = 7r/2 
for X ^ ^ while 77>(— |x|) satisfies 



tan 



^^-lg-2K(|x|) 



(12) 



according to Eq. Hll|) . For 7 7^ the phase 77> = — 3Tr/2- 
27rn at x ~ —\x\ —00. The amplitude factor for x = 
— |x| is found from Eqs. © and Hll|l 



C>(-|x|)=if(|x|) 



1 



In 



7^ + e-'^^(l^l) 



For X — |x| it is simply C>(+|2;|) — The other 



solution 



-'1:2 



(x) 



^1,2 



(— x) grows at X ^ +c 



The particle transmission Df, and hole (Andreev) re- 
flection Rfi probabilities, 1 ~ Rh + D^., are determined 
such that De — |wi(z ~ d)/wi{z = 0)|^ provided there 
are no transmitted holes, W2 = at z — d. We denote 
the x-coordinates of the end points of the trajectory at 
z = and z — d as x- and x+, respectively, such that 
dtanO = x+ — x_. For trajectories crossing the vortex 
axis, < xj_ < dtan^, x_ = — |x_|, we find 



where 



(7' + a') \osh-^[K(\x+\)+K{\x-\)] (13) 



a = cosh [K{\x+\) - if(|x_|)] /cosh [K{\x+\) + K{\x-\)] 

is the half-width of the level 7 = proportional to the es- 
cape rate of excitations from the level through the gapped 
region far from the vortex. For 7—^0, the transmission 
coefficient becomes Dp ~ cosh^^ [/f (|x+|) — iir(|x_|)]. 
It is Dp = 1 for resonant trajectories that go through 
the middle of the vortex, |x_| = jx+j = {d/2)tand. 
For trajectories that do not cross the vortex axis, 
x_,x-|_ < or x_ > 0,x_|_ > dtanO, we find Dp = 
cosh~^ [K{x+) — K{x-)] where x+ — x_ = dtan^ . 

Equation ((T^ suggests that the largest contribution to 
the transmission through the energy gapped region comes 
from the resonant tunneling when the energy coincides 
exactly with one of the levels in the vortex core. The 
transmission probability is then unity if the trajectory 
crosses the vortex close to the half of its length. How- 
ever, the width of the resonance a is exponentially small 
which leads to a small number of transmitted particles. 
This exponent, however, corresponds to only a half of the 
slab thickness, not to the entire slab thickness as it would 
be without vortices. The exponential decay of the num- 
ber of transmitted particles disappears for trajectories 
which go almost parallel to the vortex axis because the 
projection of the trajectory on the plane perpendicular 
to the vortex axis shrinks as ^ 0. The transport is thus 
determined by trajectories almost parallel to the vortex 
axis; its exponential dependence on the slab thickness is 
replaced with a power-law behavior. 
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The heat current - The energy current along z is 

dkz 



2'Km 



hkz A. 

c 



Particles u*u with the distribution n(e) carry the energy 
+e while the holes w with the distribution 1 — n{~e) 
carry the energy — e. If the electrodes are in equilibrium, 
1 — n{—e) — n{e) in each electrode. 

Consider first an infinitely thick slab. Using Eq. (O 
the energy current between the two electrodes becomes 



E 



dk: 2^ 



Excitations with positive group velocity h^^de^/ dk^ 
have the distribution ni — [e'^^^^ + as in the electrode 
1. For those with negative group velocity the distribution 
is n2 — [e*^/-^^ + I]~^ as in the electrode 2. Therefore 



Pz>0 



de 



dkz 
2'Kh 



(14) 



By the order of magnitude, the heat current is Is ^ 
{T/h){kF0[T/Tc){Ti-T2) with the thermal conductance 



Trajectories that do not cross the vortex but pass within 
an angle ~ ^l<P' near its axis give the same or- 
der of magnitude. The thermal conductance is k ^ 
{T /h){kp£)'^{£^/d)^ . The semi-classical approach requires 
krX ^ 1, which restricts the angle by » (kpS,)^^^^- 
a result, Eq. applies for d/^ ^ {kpS^Y^^. However, 
transition to the Landauer equation H15|l occurs already 
at smaller d. Thus the crossover region has be considered 
separately, taking into account non-quasiclassical correc- 
tions to the wave functions. 

In conclusion, we have shown that with respect to the 
single-electron transport the vortex core in a thick clean 
superconductor is equivalent to a mesoscopic conduct- 
ing channel where the conductance is determined by the 
Landauer formula with CdGM states playing the role of 
transverse modes. For a finite slab thickness d, the vor- 
tex core behaves like an Andreev wire: the thermal con- 
ductivity along it drops off as d~^ due to a finite order 
parameter everywhere except at the vortex axis. 
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{T/h){kpO{T/T,^ 



(15) 



in compliance with Eq. 1^)1. 

For a finite slab thickness, can be expressed through 
the transmission and reflection coefficients, 

= ivp f — \vz I / dx^ db [ Df, e (ni — 712) de , 
Jv,>o 47r 7 J 

where the sum is over all the trajectories; vp is the single- 
spin density of states at the Fermi level. The first two 
terms in the upper line are due to incoming particles and 
Andreev reflected holes on one side of the slab. The third 
term is due to transmitted particles from the other side. 

Using Eq. H13|) we find that for the resonant trajec- 
tories — dtan^ < a;_ < the integration over de- 
selects \x-\ close to The angles along the vor- 
tex core axis are small such that fc^ = kp sin 6 <^ kp 
and K{x) — {m/2h'^kr)A'{0)x^ . The localization ra- 
dius of the wave function in the x direction is A = 
S(fcF/7rmA'(0))^/^ 6(1/2 _ ^5,1/2 Let us put xq = 
— l^+l) /2 while + ~ dtanO. For smaU 
e we have K{\x+\) - -fs:(|a;_|) = xq {md/h'^kp)A' (0) 
and K{\x+\) + K{\x-\) 9 {md"^ /ti^kp)A' {Q) ~ d'^O/^^. 
Therefore, xq ^ ^'^/d and 9 ^ {^/d)^. As a result. 



h^vph\%d~'' (A'(0))~' (t2 



(16) 
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